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Abstract 

A general procedure which defines a partial ordering of cyclic permutations 
induced by continuous maps is known for constructing immediate successors 
to a cycle. We expound on this procedure in terms of labelled digraphs and 
characteristic polynomials then apply this study to period doubling, the most 
common route to chaos for a nonlinear dynamical system. 
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1 Introduction 

Li and Yorke [10 j showed that the presence of a periodic orbit with least period 3 
forces the presence of periodic orbit with all possible least periods. However before 
[10 , Sharkovsky [6, 9 had already obtained the results which include the statement 
above as a corollary. Baldwin |3] extended this idea to get a partial order on the set 
of cycles in which the type of the periodic orbit is taken into account. Bernhardt 
[5] showed how to construct immediate successors to a given cycle. 

First we reexamine [5] to obtain graph theoretical explanation of the relation 
between a cycle and its immediate successors. With this view point, we then obtain 
the representation of the period-doubling bifurcation phenomena by the digraphs 



consisting of signed vertices and directed edges and by the characteristic polynomials 
of the adjacency matrices. Obtaining a discrete expression of the period-doubling 
bifurcation phenomena implies that we may define the period-doubling bifurcation 
phenomena on the graph itself without referring to the original one-dimensional 
continuous maps. 

In section 2 we present basic definitions and concepts on the forcing relation on 
cycles in dimension one and on graph theory. In section 3 we obtain the relation 
between the digraph and the characteristic polynomial of any cycle and that of its 
immediate successor. In section 4 we describe the period-doubling phenomena of 
the family of the logistic curves with their digraphs and characteristic polynomials. 

2 Preliminaries 

Definition 2.1. Let / : / — > / be a continuous map of a compact interval I to itself. 
We define f°(x) = x and for n £ N \ {0}, f n (x) = /(/ n_1 (^))- If there exists k 6 N 
such that f k (x) = x, then we say x is periodic for / and x has (least) period s, 
where s is smallest element of N such that f s (x) = x. For x £ I, the orbit of x is 
the set {f n (x)\n > 0}. If x is periodic with least period s, then the orbit of x is the 
finite set 0{x) = {x, /(*), f(x), . . . , f^Hx)}. 

Definition 2.2. A cycle (cyclic permutation) of length n is a bijection 9 : {1,2, . . . , n} 
— ► {l,2,...,n} such that 9 k (l) ^ 1 for 1 < k < n. We denote a cycle 9 by 
9 = (ki, &2, ■ ■ • , k m ), where 9{ki) = k{ + i and 0(k m ) = k\. We assume, without loss 
of generality, that k\ = 1. Write the elements of a periodic orbit O in increasing 
order: x\ < x<i < ■ ■ ■ < x s . We say O has orbit type 9 if 9 is a cycle of length s 
and for each i G {1,2,..., s}, f(xi) = xqu\. 

(S n , o) denotes the group of permutations on n objects, C n denotes the subset of 
S n consisting of cycles (cyclic permutations) of length n. We also denote U n >i ^ n 
by C. 

Definition 2.3. The forcing relation on cycles is defined as follows: (3 £ C 
forces a £ C if and only if every continuous map of the interval that has a periodic 
orbit of type ft has a periodic orbit of type a. 

Baldwin [3] showed that this relation induces a partial order on the set of cycles 
and provided an algorithm to decide when one cycle forces another. 

Definition 2.4. [5j We say (3 £ C is an immediate successor to a £ C, or a is 

an immediate predecessor to (3, if (3 ^ a and the set of all cycles forced by (3 is 
same as the set of all cycles forced by a and a itself. 
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Definition 2.5. For 9 £ S n , the canonical ^-linear map, or the connect-the- 
dots map of 9 is defined by Lg : [1, n] — * [l,n] where Lg = 9 on {1, 2, . . . , n} and 
Lg is linear on Ii := [i, i + for i G {1, 2, . . . , n — 1}. 

It is known that [3 £ C forces a £ C if and only if has a periodic orbit of 
cycle type a [2] . 

Definition 2.6. For 9 £ C n , n > 3 we say 9 is +k-modal cycle or —k-modal 
cycle if the number of all the local extrema of Lg is equal to k and 9(1) < 9(2) or 
9(1) > 9(2), respectively. In particular +l-modal cycle is called unimodal cycle. 

The forcing relation defined by Baldwin induces a total order on the set of 
unimodal cycles [7J. 

Definition 2.7. From Lg, we define the labelled digraph (directed graph) of 

9, denoted by G(9), as follows: G(9) = (V,£) consists of the set of n — 1 signed 
vertices V = {v f, vf, . . . , v n _±} where we take vf if 9(i) < 9(i + 1) or «r otherwise, 
and the set of directed edges £ = {%•} where a directed edge from Vi to vj exists 

if Ij C Lg(Ii). 

Conversely, from G(9) we can obtain Lg. In this sense, we say they are equivalent. 
Indeed Baldwin's algorithm [3j which decides whether one cycle 9 forces another or 
not, is derived from G(9). 

Definition 2.8. The induced matrix of 9 or the adjacency matrix of G(9) [8], 
denoted by M(9) = aij , is the (n — 1) x (n — 1) matrix with (i, j)-th entry given 
by 



1, if C Lg(Ii), 

0, otherwise, 



where i,j G [l,n— 1]. 



Definition 2.9. |8 A graph H = (W, J-) is said to be a subgraph of the graph 
G = (V, £ ) if W C V and T C £. In particular if T consists of all the edges from £ 
which connect the vertices from W, then H is called an induced subgraph. 

Definition 2.10. [H] A cycle of length n in the digraph G is a subgraph with the 
vertex set {vi, . . . , v n } having directed edges from V{ to Vi+i where i = 1, . . . , n — 1 
and a directed edge from v n to V\ . In particular if v\ < Vi < ■ ■ ■ < v n then this cycle 
induces the cyclic permutation 9 = (12 ■ ■ ■ n). 
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3 Theorem 



Definition 3.1. [S] Let 9 G C n , then 0* G S2n is denned by 

0*(2A;) := 26>(fc), 0* (2k - 1) := 20(fc) - 1, 
where k € {1, 2, . . . , n}. 

Definition 3.2. [5j Let p s denote the transposition (2s — 1,2s), where s £ N. Then 
for 9 G C n define 77 as follows: 

7? := 0* op i2 o ••• op im , 

V v ' 

m 

where 1 < i\ < 12 < • ■ ■ < i m < n. 

If m is odd, then 77 G Cm and it is an immediate successor to 9. If m is even, 
then 77 G S271 \ C*2n- There are 2 n_1 immediate successors 7/ G C2 n to 0, and 2" _1 
non-cyclic permutations 7/ G S^n- (See [1] and [5]) 

For 9 and 77 above, denote their labelled digraphs by G(9) and G(tj), respectively. 
Then G{9) = (V, S) consists of the set of n— 1 signed vertices V = {v f, vf, . . . , v^^}, 
where we take tv" if < 0(i + 1) or otherwise, and the set of directed edges 
£ = {eij} where a directed edge from V{ to exists if Ij C Lg(Ii). G(rj) = (W, F) 
consists of the set of 2n — 1 signed vertices W = , wf, . ■ • , w^n-ii' w here we 
take wf if 77(7) < 77(7 + 1) or otherwise, and the set of directed edges T = {fij} 
where fij a directed edge from Wi to Wj exists if Ij C L v (Ii). Let us divide W into 
two subsets Wi = {wf, • • • , ^2^-2} an< ^ ^2 = { w f, wf, . . . , u^n-l}* 

Then we have the following observations: 

(7) G(t7)|wi, the induced subgraph of G(rj) = (W,.F) determined on Wi is iso- 
morphic to the digraph G{9) = (V,£). Indeed, the isomorphism (p : V — ► Wi is 
given by 

= W 2i 

where i G {1, 2, ... ,n — 1}. Observe that the sign of vertex Tjj is always consistent 
with that of vertex w<ii- 

{%%) In G(?7) = (W, .F) there does not exist any directed edge from W2 to W\. 

(in) G(r])\w 2 , the induced subgraph of G(rj) = (W,F) determined on W2, con- 
sists of an Ti-cycle induced by 9. Indeed, this n-cycle is given by 
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^20(1)-! ^20 2 (1)-1 



^20"- i(l)-l 



and the type of this n-cycle is also 0. 

Together, (i), (ii), and (Hi) provide a graph theoretical explanation of the reason 
why r\ is an immediate successor to 9 (see Lemma 2.4 in [5]). As a consequence, we 
have the following theorem. 

Theorem 3.1. Denote the characteristic polynomials of the adjacency matrices of 
G(9),G(r]),G(r))\y\; 1 and G(r])\w 2 by Pg(X),P Tt (X),Pw 1 (X) andP\\> 2 (X), respectively. 
Then P V (X) is given by the product of Py\; 1 (X) and Pyy 2 (A) ; namely 



As a corollary of this theorem, we have the following. 

Corollary 3.1. The spectrum of G(rf) consists of the spectrum of G(0) and the n 
distinct roots of unity. 

Assume that the adjacency matrix of G(6) is given by 



P V (X) = Pg(X) ■ (X n - 1). 



M(6) 



- (n-l)x(n-l) 



where 1 < 



i, j < n — 1) and the adjacency matrix of G(rj) is given by 



M(rj) 




. (2n-l)x(2n-l) 



where 1 < i,j < 2n — 1. 

Then we have the following Lemma 3.1, and 3.2. 



Lemma 3.1. G(?7)|wi is isomorphic to G(6) by <p, namely 



b2i,2j — a>i,j, 



where 1 < i,j < n — 1. 
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Proof. For any 9 £ C n , 9* <E S 2n is defined by 



2i - 1 2i 2i + l 2i + 2 

20(0-1 29{i) 2(9(2 + 1) - 1 20(i + 1) 



If 77 := 9* o ■ ■ ■ o pi o ■ ■ ■ , then 
rj = 



2i-l 2i 2i + l 2i + 2 

20(0 20(0-1 20(i + l)-l 20(i+l) 



and 

[20(0 - 1, 26»(i + 1) - 1] if 6{i) < 9{i + 1) 

[20(i + 1) - 1, 20(i) - 1] if 0(i)>0(i + l) ' 

If 77 := 9* o • • • o p i+1 o ■ ■ • , then 



L„([2i,2t + 1]) 



77 

and 



2«-l 2i 2i + l 2i + 2 

2(9(z) - 1 20(«) 26>(i + 1) 2(9(i + 1) - 1 



r rr%- 9,- + m- J [2<?(i),20(* + i)] if 0(i)<0(* + i) 

1,[[ l, l+ l) -\[29(i + l),29(i)} if 0(i)>0(i+l) 
If 77 := 9* o • • • o pi o o • • • , then 



77 

and 



2i-l 2« 2« + 1 2i + 2 

20(0 26»(i) - 1 29{i + l) 20(i + 1) - 1 



r no; 9j + m-i [20(i)-l,20(i+l)] if 0(0 < 0(i + 1) 
" U ' + Jj " \ [20(i+l),20(i) -1] if > + 1) 

If r) := 9* o p ix o p i2 o ■ ■ ■ o p im does not include nor pi+i, then 



77 

and 



2i-l 2i 2i + l 2i + 2 

2(9(«) - 1 29{i) 20(t + l)-l 29{i+l) 



r n 0j o,- + in_/ [20(i),20(i + l)-l] if 0(i)<0(* + l) 
T ' l[ ' + Jj ~ \ [2©(i+ 1)- 1,20(0] if 0(*)>0(» + l) 

Under the assumption 9{i) < 9(i + 1), we have 

a-ij = 1 
& Ij C L fl (/i) 

[j,j + l]c[9(i),9(i + l)} 
& 9{i) < j and j + 1 < 9(i + 1) 
4^ 20(i) < 2j and 2j + 1 < 2(9(i + 1) - 1 
& [2j,2j + l] C [20(0, 20(7; + !)-!]. 
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Here we note that 



[20(i),20(i + l)-l] c [20(t)-l,20(i+l)-l],[20(t),20(i+l)] c [20(i)-l,20(i + l)], 
and 



[ 2j , 2j + 1] ± [29(f) - 1,26(1)}, [29(i + 1) - 1, 29(i + 1)], 

even odd odd even odd even 



thus, 



[2j,2j + 1] C [20(i),20(i + l)-l] 
& [2j,2j + l] C [r l (2i),r ] (2i + l)} 

44> l2j C L v (l2i) 
^ &2i,2j = 1- 

Similarly, under the assumption 9(i) > 9(i + 1), we can also show that 



(Hj — 1 &2i,2j 

This concludes the proof of the lemma. 



1. 



□ 



Lemma 3.2. In G{rj) = (W, J 7 ), there does not exist any directed edge from W2 to 
Wi, i/iai is 

°2i-l,2j = 0, 

where 1 < i < n, 1 < j < n — 1. Moreover, G(n)\w 2 consists of an n-cycle induced 
by 9, namely 



&2i-i,2j-i = 1 ifj = 0{i), 
&2i-i,2j-i = otherwise, 



where 1 < i,j < n. 
In other words 



b\2 

&32 



&14 
^34 



&(2ra-l)2 &(2n-l)4 



&l(2n-2) 
&3(2n-2) 



&(2ra-l)(2n-2) 



nx(n-l)i 



and 



611 


613 


b\{2n-\) 




" Cn 


C12 • 


■ Ci n 


&31 


&33 


&3(2ra-l) 




C21 


C22 • 


■ c 2n 


&(2n-l)l 


&(2n-l)3 • • 


• &(2n-l)(2ra-l) . 






Cn2 ■ 


C-nn 
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&2i-l,2j-l — Cij, 

where 

Cij = l ifj = 0(i), 
Cij = otherwise. 

Proof. 

3 = 0(i) 

^ 6*(2i - 1) = 2j - 1 and 9*{2i) = 2j 
^ [2j-l,2j] = [0*(2i-l),0*(2i)]. 



2^ — 1 2z 
*(2i-l) 0*(2i) 



For any 6> G C n , 0* G S2 n is denned by 
0* := 

If 77 := 0* o • • • o pi o ■ ■ ■ , then 

_ ( ■ ■ 2* — 1 2z 
r/ ~ V ••• 0*(2z) 0*(2i-l) 

and thus [2j - 1, 2j] = [9*(2i - 1), 6*(2i)} = [77(2*), fj(2t - 1)]. 
If r/ := 0* o p ix o pj 2 o • • • o p im does not include pi, then 



7] 



2i-l 2i 
0*(2i-l) 9*(2i) 



and thus [2j - 1, 2j] = [0*(2i - 1), 9* (2i)} = [q{2i - l),»/(2i)]. 
In either event, we have 

[2j-l,2j] = [0*(2z-i),0*(2*)] 

^ ^2j-l = L v (l2i-l) 
b2i-l,2j-l = 1 , 

b 2i -i, k = if A; ^ 2j - 1, 



44> 



ending the proof of the lemma. □ 
The following statement is obvious: 



Lemma 3.3. 

Cmxl B mxm 

We are now ready to prove the theorem. 



\A lxl \\B. 



mxm\ ■ 



8 



Proof of Theorem 3.1. The characteristic polynomial of the adjacency matrix of 
G{9) is given by 



Pe(X) = \\I-M(0)\ 



A - an -a 12 
— a2i A — a22 



— a (n-l)l _a (n-l)2 



A- 



Let 



C2n-1 : = 



1 2 3 



n n + 1 n + 2 n + 3 



-«l(n-l) 
-«2(n-l) 



a (n-l)(n-l) 



In- 1 



1 3 5 ■■■ 2n-l 2 4 6 ■ ■ ■ 2n - 2 
Using Lemmas 3.1, 3.2, and 3.3, the characteristic polynomial P V (X) of n is thus 



X-b n 
-hi 



-bi2 
X - b 2 2 



-°(2n-2)l —b(2n-2)2 
-b(2n-l)l —b(2n-l)2 



—h(2n-2) 
—b2(2n-2) 



-6 



? l(2n-l) 
'2(2n-l) 



A-6 



(2n-2)(2n-2) 



?(2n-2)(2n-l) 



Sgn(cr 2n _i) 



A-6 



ii 



>12 



-6, 



(2n-2)l — °(2n-2)2 



^21 



A-6. 



22 



(sgn(cr 2n -i)) 5 



b(2n-l)l —b(2n-l)2 

X-bn 



b(2n-l)(2n-2) X — 6(2n-l)(2n-l) 

— &l(2ra-2) _ &l(2n-l) 



A-6 



(2n-2)(2n-2) 
l) 2{2n-2) 



5 (2n-l)(2n-2) 



5 (2n-2)(2n-l) 



5 2(2n-l) 



A-6 



(2n-l)(2n-l) 



"°(2n-l)l 
-621 

-&(2n-2)l 



-6 



l(2n-l) 



-6l2 



A — b(o n _ 



(2n-l)(2n-l) ~~ &(2ra-l)2 



A -611 



-6 



l(2n-l) 



-6 



(2n-l)l 



X — 6( 2n -l)(2n-l) 



- & 2(2n-l) 

&(2ra-2)(2n-l) 
A-6 2 

-6r 



r ^22 



^(271-2)2 



A - 6 2 2 

~b(2n-2)2 

-b 

■■ A-6 



-6 



l(2n-2) 



-fr(2n-l)(2n-2) 
— ^2(2n-2) 

A _ &(2n-2)(2ra-2) 

'2(2n-2) 



(2n-2)(2n-2) 



(X n + sgn(0)(-l) n )P 9 (X) 
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since (sgn(a 2 n-i)) 2 = 1, 



"0(2n-l)2 ' ' ' — 0(2n-l)(2n-2) 



rax (n— 1) i 



A - 6 2 2 

~&(2n-2)2 



5 2(2n-2) 



A - 6, 



(2n-2)(2n-2) 



and 



A-6 
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A - an 

— a (n-l)l 
' ? l(2n-l) 



-ai( n _i) 



A — a 



(n-l)(n-l) 



— &(2n-l)l " ' ' A — 6(2n-l)(2n-l) 

A - cn • • • -ci„ 

C-nl ' A C nn 



where 



Cij = 1 if J = ^(i) 
= otherwise. 



sgn(t) • (A-O) n + sgn(0) • (-1) 
A n + sgn(#)(-l) n . 



Here we note that sgn(#)(— l) n is always equal to —1 since 

sgn(0)(-l) n = 
This completes the proof of the theorem. 



(+l)-(-l) if n is odd, 
(— 1) • (+1) if n is even. 



4 Immediate Successors and Period Doubling 

4.1 The logistic map 

The logistic map [TJ [9] is written as 



□ 



(l) 



where the current value x n , is mapped onto the next value x n+ \. Here, we restrict 
our considerations to < x n < 1 and < a < 4. The corresponding function to the 
logistic map, called the logistic function or the logistic curve, is given by 



f a (x) = ax(l - x) 



(2) 
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By repeatedly iterating the logistic map forward through time, we may observe 
different behaviors of the iterated solutions. The sequence of iterated solutions to 
the map is called an orbit. The behavior of orbits, which originate from typically 
chosen initial condition used for the iterations, depends on the control parameter a. 
The logistic curve is parabolic like the quadratic function with / o (0) = / a (l) = 
and a maximum at x = 1/2. Being a smooth curve with only one critical point, 
that is, a local maximum, it is classified as a unimodal (single humped) map. 



4.2 Period-doubling bifurcations 

The logistic map has two fixed points: and {a — 1) / a. If a < 1, then 
I fa (0)1 < 1 an d we say that the fixed point x = is stable for a < 1. If a > 1, then 

1 fa (0)1 > 1 an d so the fixed point x = is now unstable and a new stable fixed point 
x = (a — l)/a emerges. The fixed point x = (a — l)/a is stable for 1 < a < 3. At 
a = 3, the fixed point x = (a — l)/a becomes unstable and a stable 2-cycle (period 

2 attractor) {x,y} = {(a + 1 — V a 2 — 2a — 3)/2o, (a + 1 + vV - 2a — 3) /2a} is 
said to bifurcate from this fixed point. The above-mentioned 2-cycle is the solution 
of the system of equations given by y = ax(l — x) and x = ay(l — y). When 
a = 1 + v6 ~ 3.44949, this 2-cycle loses its stability and a stable 4-cycle begins 
followed by a stable 8-cycle (between a ~ 3.54409 and a ~ 3.56440). The period 
continues doubling over shorter intervals until a ~ 3.56995 where the chaotic regime 
takes over. A period 3 attractor is born at a = 1 + y/8 ~ 3.82843 [11], causing a 
notably large window in the bifurcation diagram (see Figure 1). This 3-cycle also 
undergoes a period doubling cascade in which period 3-2^ attractors are successively 
produced until chaos. This mechanism by which each (fc2^)-cycle loses its stability 
and is replaced by a stable (fc2^ +1 )-cycle is called period-doubling bifurcation 
(see Figure 2). 




i -i 



Figure 1. Bifurcation diagram of the logistic map 
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0,3 3.54 3.545 3.55 3.555 3.56 3.565 3.5T 



Figure 2. Magnification of the first few period doubling 

To determine period 4, period 8, and in general, period n solutions of the logistic 
map, we solve the system of difference equations, x\ = X2, X2 = x%, . . . , and 
We note that these computations become very complicated even for relatively small 
values of n. Finding the location of these bifurcation points on the bifurcation 
diagram also becomes very difficult. Even so, examining the types of orbits and 
their characteristic polynomials turns out to be less tortuous. This is illustrated in 
the following sections. 

4.3 Period-doubling bifurcation beginning with nontrivial fixed point 

The attracting period 2 orbit of type 9\ = (12) bifurcates from the fixed point of 
type $o = (1) a t a = 3. 6\ is the only one immediate successor to 6o, and it follows 
that 

P 9l (X) = A-l 

= (A-l).P eo (A), 

if we define P$ Q {\) := 1 for 6>o = (1). Indeed the graph G{9\) = (Vi,£i) consists 
of a vertex v\\ with the minus sign and a directed edge from v\\ to itself (a loop 
associated with vu): 




The attracting period 4 orbit of type 02 = (1324) bifurcates from the period 2 
orbit of type 6i at a = 1 + V6 ^ 3.44949. 6 2 with modality +1 is one of 2 2 " 1 = 2 
immediate successors to 6\. Here we note that (1423) with modality —1 is another 
immediate successor to 9±, however, the logistic map can realize only +l-modal 
(unimodal) cycles. Now it follows that 

Pe 2 (X) = A 3 -A 2 -A + l 
= (A 2 -1)(A-1) 
= (A 2 -l)-P ei (A). 
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Indeed the graph G(02) = ^2,^2) is given as follows: 




where V2 = {^21' u 22' ^23}' ^21 = {^22}' anc ^ ^22 = {^21' ^23}- Here we note that 
(i) G(02)\\> 21 i s isomorphic to G(6i), the sign of vertex v\\ is consistent with that of 
vertex W22, (ii) there are no directed edges from V22 to V21, (Hi) G(#2)|v 22 consists 
of a 2-cycle induced by 6\ . 

The attracting period 8 orbit of type #3 = (15472638) bifurcates from the period 
4 orbit of type 62 at a « 3.54409. #3 with modality +1 is one of 2 4_1 = 8 immediate 
successors to 62- There are other 7 immediate successors to 62, however, their 
modalities are not equal to +1. Now it follows that 

p d3 (X) = a 7 -A 6 -A 5 + A 4 -A 3 + A + A-1 
= (A 4 -1)(A 2 -1)(A-1) 
= (A 4 -1)-P, 2 (A). 

Indeed the graph G(6s) = ^,£3) is given as follows: 




where V 3 = {v^,v^, v 33 , v 34 , v 35 , v 36 , v 37 }, V31 = {v£ 2 , v 34 , v 36 }, and V32 = {wjji, W33, v 35 
Here we note that (i) G(#3)|v 31 is isomorphic to £(6*2), the sign of vertex v 3 ( 2 j) is 
consistent with that of vertex V2j (j = 1, 2, 3), (ii) there are no directed edges from 
V32 to V31, (Hi) G(^3)|v 32 consists of a 4-cycle induced by 62- 
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Inductively, for the attracting period 2 orbit of type 0£ with modality +1, which 
is one of 2^ 2 -1 ) immediate successors to 0£-±, and bifurcated from the period 2^ _1 
orbit of type 0£-i, we obtain G{0£) and Pe e (X). 

4.4 Period-doubling bifurcation beginning with period 3 

Next consider the period-doubling bifurcation beginning with period 3 orbit of 
type 0o = (123). The graph G(0q) = (Vo,£o) consists of two vertices with a 2-cycle 
and a loop: 




and thus Pe (X) = A 2 - A - 1. 

The attracting period 6 orbit of type 0\ = (135246) bifurcates from the period 3 
orbit of type #0 = (123) at a « 3.8415. 0\ with modality +1 is one of 2 3 " 1 = 4 im- 
mediate successors to 0q. (146235), (136245), (145236) are also immediate successors 
to 0o, however, their modalities are not equal to +1. Now it follows that 

p 9i (\) = a 5 -A 4 -A 3 -A 2 + A + 1 
= (A 3 -1)(A 2 -A-1) 
= (A 3 -l).P eo (A) 

Indeed the graph G(0\) = (Vi,£i) is given as follows: 




where Vi = Ki,^,^^^,^}, Vn = {vf 2 ,v u }, and V 12 = K,^ 3 ,%}. Here 
we note that (i) G(^i)|v n is isomorphic to G(0q), the sign of vertex t>i(2j) is consistent 
with that of vertex voj (j = 1,2), (ii) there are no directed edges from V12 to Vn, 
(in) G(0i)\v 12 consists of a 3-cycle induced by ^o- 
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Inductively, for the attracting period (3 • 2 ) orbit of type 0£ with modality +1, 
which is one of 2^' 2t immediate successors to and bifurcated from the 

period 2 e ~ l orbit of type we obtain G(9() and Pg e (X). 

4.5 Period-doubling bifurcation beginning with period k 

It is known that periodic attractors exist for all odd integers k. Therefore, in 
general, for the attracting period (k ■ 2 £ ) orbit of type 0£ with modality +1, which is 
one of 2^ k ' 2 immediate successors to Op-i, and bifurcated from the period 2 e ~ 1 
orbit of type 0£—i, it follows that 

£ 

1=1 

where Pg (X) is the characteristic polynomial of the induced matrix of 9$, and 9q is 
the type of the original period k orbit. 

Indeed the graph G{0^) consists of the set of vertices Vi and the set of directed 
edges S e where Ve = {«a, • , v l{k . 2 i_ 2) , v l{k . 2l _ x) \, Vtl = {vn, ■ ■ ■ ,v e{k . 2 i_ 2) }, 

and V12 = {fii,%,-<M^(MM)}- Here we note tnat (*) G(6i)\y n = G($i-i) for 
i = 1,2, ... ,1, and the signs of the vertices are preserved by each of the above iso- 
morphisms, (ii) there are no directed edges from V-a to Vn for i = 1, 2, . . . ,£, (in) 
G(6i)\y i2 consists of a (k ■ 2 l ~ 1 )-cycle induced by 6i-i for i = 1, 2, . . . , £, and G(0q) 
consists of k — 1 vertices and some directed edges. 

In other words, G(9e) essentially consists of k — 1 vertices with some directed 
edges, a fc-cycle of type Qq, a 2£;-cycle of type 9±, ■ ■ ■ , a (k ■ 2^ -1 )-cycle of type 
That is the reason why P$ t (x) can be expressed as a product of Pg (X), X k — 1, 
X 2k — 1, • • • , X^ k ' 21 1 - ) — 1. Thus the spectrum of G(6e) consists of the spectrum of 
G(9 Q ) and (k ■ 2^) distinct roots of A^' 2 * -1 ) -1 = where i = l,2,...,£. 

Remark 4.1. Conversely, we may define G(6t), the doubling of the graph G(0j_i) 
by properties (i), (ii), and (Hi). 

5 Acknowledgement 

The authors wish to thank Professor Evelyn L. Tan for referring to us [8]. 

References 

[1] Paul S. Addison, Fractals and Chaos An Illustrated Course, IOP Publishing Ltd 
(1997). 



15 



[2] L. Alseda, J. Llibre, M. Misiurewicz, Combinatorial Dynamics and Entropy in 
Dimension One, 2nd ed., World Sci., (2000). 

[3] S. Baldwin, Generalizations of a theorem of Sarkovskii on orbits of continuous 
real-valued functions, Discrete Mathematics, Vol 67 Issue 2 (1987), 111-127. 

[4] C. Bernhardt, Simple permutations with order a power of two, Ergod Th. & 
Dynam. Sys., 4, 179 - 186 (1984). 

[5] C. Bernhardt, The ordering on permutations induced by continuous maps of the 
real line, Ergod Th. & Dynam. Sys., 7, 155 - 160 (1987). 

[6] K. Burns, B. Hasselblatt, Sharkovsky's Theorem, (preprint), (2007). 

[7] P. Collet and J.-P. Eckmann, Iterated Maps on the Interval as Dynamical Sys- 
tems, Progress in Physics, vol. 1, Birkhauser, Boston, Mass., (1980) 

[8] D. M. Cvetkovic, M. Doob and H. Sachs, Spectra of Graphs: Theory and Appli- 
cations, (3rd rev. enl. ed.), New York, Wiley, (1998). 

[9] P. G. Drazin, Nonlinear systems (Cambridge Texts in Applied Mathematics), 
Cambridge University Press (1997). 

[10] T.-Y. Li, J. A. Yorke, Period three implies Chaos, Amer. Math. Monthely, 82 
no.10 (1975), 985-992. 

[11] P. Saha, S.H. Strogaz, The Birth of Period Three, Mathematics Magazine, Vol. 
68 No 1 Feb. (1995), 42-47. 



16 



